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ABSTRACT 


The maximum-flow network-interdiction problem (MXFI) arises when an 
interdictor, using limited interdiction resources, wishes to restrict an adversary’s use of a 
capacitated network. 

MXFI is not easy to solve when converted to a binary integer program. Derbes 
(1997) uses Lagrangian relaxation to solve the problem, at least approximately, for a single 
value of available resource, R. Bingol (2001) extends this technique to solve MXFI 
approximately for all integer values of 7? in a specified range. But, “problematic 7?-values” 
with substantial optimality gaps do arise. We reduce optimality gaps in two ways. First, 
we find the best Lagrangian multiplier for problematic 7?-values by following the slope of 
the Lagrangian function. Secondly, we apply a limited-enumeration branch-and-bound 
algorithm. 

We test our algorithms on six different test networks with up to 402 nodes and 1826 
arcs. The algorithms are coded in Java 1.3 and run on a 533 MHz Pentium III computer. 
The first technique takes at most 39.3 seconds for any problem, and for one instance, it 
solves 8 of the problem’s 15 problematic 7?-values. For that problem, the second technique 
solves four of the remaining problematic 7?-values, but run time increases by two orders of 
magnitude. 
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EXECUTIVE SUMMARY 


This thesis is concerned with solving or approximately solving a maximum-flow 
network-interdiction problem (MXFI). This problem arises when an interdictor, using 
limited interdiction resources, wishes to restrict an adversary’s use of a capacitated 
network. In particular, the adversary, or “network user” seeks to maximize flow through 
the network, and the interdictor seeks to minimize that maximum flow by destroying arcs. 

MXFI is not easy to solve when converted to a binary integer program. Derbes 
(1997) uses Lagrangian relaxation to attempt to solve this problem for a single value of 
interdiction resource, R. Bingol (2001) extends this technique to solve MXFI for all integer 
values of R in some specified range assuming that exactly one unit of resource is required 
to interdict an arc. However, about 25% of these values are “problematic” in tests, i.e., 
positive optimality gaps arise. 

In our study, we first focus on decreasing the optimality gaps by finding the best 
Lagrangian multiplier. We simply follow the slope of the Lagrangian function to 
accomplish this. For most test problems, this method decreases the gap for problematic R- 
values by approximately 50%. This method helps find better feasible solutions and 
sometimes prove optimality. But, some large optimality gaps may remain. 

We tackle the remaining problematic 7?-values with a restricted branch-and-bound 
technique, which has only one level of enumeration below the root. This technique creates 
a strict partition of the feasible region and the minimum lower bound from among these 
partitions is a lower bound on the optimal solution value. In our tests, we solve 
approximately 70% of the problematic 7?-values that remain after applying the first 
technique. However, branch-and-bound increases run times dramatically. 

Both procedures are written and compiled using the Java 1.3 programming 
language. All tests are performed on a personal computer with a 533 MHz Pentium III 
processor and 192 MB of RAM, running under the Windows XP operating system. For test 
purposes, we have used 6 different grid networks. The smallest one has 27 nodes and 86 
arcs and the largest one has 402 nodes and 1826 arcs. 
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Applying both techniques together, we have been able to either solve problematic 
i?-values exactly or improve their optimality gaps for most cases. In fact, 85% of the 
problematic i?-values are solved optimally in our tests. One problem with a large 
optimality gap does not improve, however. 



I. INTRODUCTION 


This thesis develops a method based on Lagrangian relaxation to solve a maximum- 
flow network-interdiction problem. This problem arises when an interdictor, using limited 
interdiction resources, wishes to restrict an adversary’s use of a capacitated network. In 
particular, the adversary, or “network user” seeks to maximize flow through the network, 
and the interdictor seeks to minimize that maximum flow by destroying arcs. 

The maximum-flow network-interdiction problem (MXFI), converted to an integer 
program, is difficult to solve because of an interdiction-resource constraint. Derbes (1997) 
uses Lagrangian relaxation to move this constraint into the objective function, and is able 
to solve MXFI efficiently for certain resource levels R. However, he solves for only one 
value of 7? at a time. The goal of this thesis is to solve MXFI efficiently for all reasonable 
values of 7? in a single procedure. 

The underlying reasons that motivate us to obtain solutions for all values of R are; 
First, R may not be specified. For instance, at an early stage of planning, we may not know 
how much resource will be available when the interdiction plan is to be carried out. Thus, 
we would like to provide the decision-maker with solutions for all potential values of R. 
Secondly, decision-makers often like to have some input on the amount of resource to be 
expended. This helps them review and maybe revise their decisions. Finally, it is not 
much harder to solve MXFI for all reasonable values of R than for a single value. 

Derbes models the interdiction of arc k as requiring a small, integer amount of 
resource, > 0. For the sake of simplicity, we assume = 1 for all interdictable arcs k, 

although arcs may also be specified as non-interdictable. Derbes also considers a dynamic 
version of MXFI in which flows require time to move through the network and interdicted 
arcs can be repaired over time. The dynamic version of MXFI is beyond the scope of this 
thesis. 

Bingol (2001) investigates our version of MXFI and also attempts to solve for all R 
using a Lagrangian-relaxation procedure. Essentially, he estimates the values of the 
Lagrangian functions for all R simultaneously. He assumes that these functions only break 
at points at which the Lagrangian multiplier equals the capacity of some arc. In carrying 
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out computations, he encounters some values of R for which his procedure cannot find an 
exact solution, so he develops a heuristie to handle these cases and determines assoeiated 
optimality gaps. He does not, however, optimize the Lagrangian-based lower bounds. 

In this thesis, we use a proeedure similar to Bingol’s to find solutions for all values 
of R but will optimize the Lagrangian lower bounds as needed. When signifieant 
optimality gaps are identified, we apply a restrieted branch-and-bound procedure to 
improve the lower bounds and possibly identify better solutions. This proeedure is 
restricted to ereate enumeration trees having limited depth. 

A, BACKGROUND 

There are different types of network-interdietion problems. For instance, Israeli 
(1999) focuses on Maximizing the Shortest Path (MXSP) to slow a network user’s 
movement between two speeified nodes in a network. The k-most-vital-arcs problem is a 
special case of MXSP in whieh the interdictor seeks to destroy exactly k arcs to maximize 
the shortest path length. Steinrauf (1991) designs a model to isolate a targeted demand 
node in a network. Morton (2001) describes a stoehastie network-interdiction model to 
help detect the smuggling of stolen nuclear materials. However, in this thesis, we are only 
eoncerned with the deterministic maximum-flow network-interdiction problem, MXFI. 

The seientific literature on network interdiction begins with Wollmer (1966) who 
foeuses on “Removing Ares from a Network,” where the objeetive is to maximize the 
reduction of the flow between an origin and a destination node. Later, Wollmer (1970) 
presents two algorithms for targeting strikes in an LOC (lines-of-communication) network 
attempting to make the network user’s costs as large as possible over time. The LOC 
networks are deseribed by a direeted network, with arcs representing road, rail, or 
waterway segments. The eosts in his study inelude both flow costs in dollars, vehicle 
hours, manpower units, or any other appropriate measures; plus repair eosts ineurred by 
strikes, measured in similar units. 

MeMasters and Mustin (1970) devise an algorithm to determine the optimum 
interdiction plan for minimizing network flow capacity when the minimum capaeity on an 
are is positive and the eost of interdiction is linear in the amount of arc capacity reduction. 
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Ratiff, Sicilia and Lubore (1975) show how to find a set of arcs, called “the N most 
vital links,” whose simultaneous removal from a single-commodity flow network results in 
the greatest decrease in the throughput capacity of the remaining system. 

Steinrauf (1991) develops two mathematical programs to determine strategies to 
interdict a network using limited resources. The first model identifies a set of arcs whose 
interdiction minimizes the maximum flow through the network, constrained by the 
available resources. This is essentially MXFI. The second model identifies a set of arcs to 
interdict that isolates a target node and the largest possible set of contiguous nodes. The 
latter model can be used if the exact location of a target node is uncertain. 

Wood (1993) shows that the MXFI is NP-complete even if the interdiction of an arc 
requires exactly one unit of resource. New, flexible integer-programming models are 
developed for the problem and a number of variants (partial arc interdiction, multiple 
sources and sinks, undirected networks, multiple resources, multiple commodities). Valid 
inequalities and a reformulation are derived to tighten the LP relaxations of some of these 
models. 

Cormican (1995) solves MXFI using Bender’s decomposition, which takes 
advantage of easy-to-solve network-flow subproblems. Cormican et al. (1998) solve 
stochastic versions of MXFI where arc capacities and/or interdiction success may be 
uncertain. Stochastic versions of MXFI are beyond the scope of this thesis, however. 

Whiteman (1999) discusses techniques for planning the interdiction of a complex 
infrastructure network to improve single-strike effectiveness. The motivation for his work 
arises from the huge cost of modern weapons and delivery platforms and intolerance for 
casualties, and the consequential desire to minimize the number of subsequent attacks. To 
optimize target selection, he modifies the integer program developed by Wood (1993), 
specifying different flow goals. 

Bingol (2001) extends the earlier work by Derbes (1997) on a Lagrangian- 
relaxation technique to solve MXFI for all integer values of total interdiction resource, R, 
in some specified range. His basic procedure solves MXFI exactly for most values of R, 
but “problematic values” of R arise with positive optimality gaps. 
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B, 


OUTLINE OF THESIS 


This chapter has provided baekground on network interdietion problems and has 
surveyed previous studies in this field. The following chapter explains MXFI in detail, and 
develops the Lagrangian-relaxation teehnique for solving it, at least approximately. In 
Chapter III, we first study the algorithms used by Derbes (1997) and Bingol (2001). 
Subsequently, we present two new proeedures to improve optimality gaps for problematie 
7?-values and sometimes even solve these problems exaetly. Chapter IV provides the 
eomputational results for the new proeedures. Finally, Chapter V provides eonelusions and 
points out directions for further researeh. 
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II. FUNDAMENTALS OF NETWORK INTERDICTION 


This chapter first defines neeessary eoneepts and notation for our study. The 
notation is similar to that used by Wood (1993) and Bingol (2001). The seeond part of the 
ehapter starts with MXFI, whieh is a min-max model, and develops a minimizing integer¬ 
programming formulation (MXFI-IP) and a Lagrangian relaxation (MXFI-LR) of that 
integer program. 

A, DEFINITIONS AND NOTATION 

G = {N, A) denotes a direeted network with node set N and are set A. An are k = {i, 
j) starts from “tail node” i and ends at “head node’’^. Eaeh are k has a capacity Uk > 0. The 
network has a souree node s and a sink node t. For mathematical purposes, we add to A the 
artificial arc a = {t, s) with u^=cc . The forward star of node i, FS{i), includes all arcs of 
the form The reverse star of node i, RS{i), represents all arcs of the form (j, i). 

In this study, only arc interdiction is allowed. In the test networks, is the amount 
of interdiction resource neeessary to destroy are k. For simplieity, we assume = I or rk 
= oo; the latter case simply means the are cannot be interdicted. 

A cut {Ns, N) partitions the node set N into two subsets. Ns and Nt where s e Ns and 
t e Nt. An are {i, j) is a forward arc of out {Ns, N) if i e Ns and j e Nt. 

Ac represents the set of forward ares for out {Ns, Nf The capacity of a cut, u{Ns, 
Nt), is the sum of the capacities of the forward ares in that cut: 

u{Ns,Nt)= • 

keA^ 

If u{Ns, Nt) is the minimum of all cut capacities, then {N^,Nf is a minimum cut. 
We know by the maximum-flow minimum-out theorem that the capacity of the minimum 
out equals the maximum s-t flow in that network; the maximum flow model will be 
described shortly. 
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B, 


NETWORK INTERDICTION MODEL 


Here, we define the baste MXFI, whieh is then eonverted into MXFI-IP. 
Lagrangian relaxation for MXFI-IP is presented subsequently. 

1, Maximum Flow Network Interdiction Problem (MXFI) 

In MXFI, an interdietor, using limited interdietion resourees, destroys ares of the 
enemy’s network G= {N, A) and makes them useless. The objeetive is to minimize the 
maximum amount of flow that the enemy, i.e., the network user, ean transport through the 
network. We define MXFI as a funetion of total interdietion resouree R: 


MXFI(i?): 

Indices: 

i,jeN 
s e N 
t e N 

k = {i,j) eA 
a = {t, s) eA 

Data: 

Uk 

ric 

R 

Decision Variables: 

yk 

Xk 


nodes of network G={N,A) 
souree node 
sink node 
network ares 

artifieial return are ineluded in A 

nominal, uninterdieted eapaeity of are k 

amount of resouree neeessary to interdiet are k {ra = go) 

total amount of available resouree 

amount of flow on are k (network user) 

1 if interdietion of are k is attempted (interdietor) 

0 otherwise 
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Formulation: 


z{R) = min max 

xeX y 


S-f-i - Hyk =0 

V/ e N 

(1) 

k&FS{i) k^RS{i) 




\/k e A 

(2) 

where X = i x e {0,1}'"^' 

Y,hx,<R,x^=o\. 

(3) 

1 

k&A j 



In this model, constraints (2) force capacities down to zero for all interdicted arcs. 
For fixed x, the inner maximizing problem is a standard maximum-flow model, which 
maximizes the flow from a source node to a sink node subject to flow-balance constraints 
(1) and arc capacities (2). MXFI(i?) can also be reformulated as MXFIR(i?) (Cormican et 
al. 1998); 


MXFIR(R): 


z{R) = mmmaxy^-^x,y, 

xeX y 


keA 


keFS{i) k^RSii) 


yi eN 


0 ^ ^ Wj 


yk e A 


where X = j x e |0,1}'“^' ^ < i?, = 0 


keA 


2 . An Integer Program for MXFI (MXFI-IP) 

Taking the dual of the inner maximization of MXFIR gives us a linear minimization 
model (which is derived by Wood 1993 directly from MXFI); 
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MXFI-IP(i?) 
Indices: As in MXFI 


Data: As in MXFI 

Decision Variables: 

tti 1 if / e Nt and 0 if / e Ns, for a model-identified eut {Ns,Nt) 

fik if A: is a forward are of out {Ns,Nt) and not interdioted, 0 otherwise 

Xk 1 if interdiotion of are k is interdioted, 0 otherwise (interdiotor) 

Formulation: 

z{R) = mmY,uJ, 


s.t. 

a.-aj+x,+p, >0 

\/k = (i,j) e A -{a} 

(4) 


+Xa +Pa 

for a = {t, s) 

(5) 




(6) 


keA 


ai e {0,1} Vi e N 
Xk, Pk^ {0,1} yk ^ A 

as = Q, at = \ 

Pa — = 0 

3. Lagrangian Relaxation for MXFI (MXFI-LR) 

Lagrangian relaxation is a general solution strategy for solving mathematioal 
programs that permits us to deoompose problems to exploit their struoture. This solution 
approaoh leads to bounds on the optimal objeotive funotion value and, frequently, to good, 
though not neeessarily optimal solutions (Abuja, Magnanti and Orlin 1993, p. 601). 
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MXFI-IP can be difficult to solve because of constraint (6), so we transfer this 
constraint into the objection function following Derbes (1997). After this modification, the 
new model, MXFI-LR, can be easily solved, although it may not always yield an optimal 
solution. The model is: 


MXFI-LR 

z{A, R) = min ^ - R) 

keA k^A 

s.t. a. -Uj + x,^ + >0 \/k = (ij) e A - {a} 


a, -a^ + x^ + P^>1 for a = (t, s) 
at G {0,1} Vi G N 
Xk, Pk ^ {0,1} yk e A 

as = 0, at = \ 

Pa = Xa = 0 


(V) 


The term Ai^Vj^x^. - R) represents a penalty for exceeding or a reward for not consuming 

keA 

the resource budget R. The parameter A can be interpreted as the dual cost of interdicting 
an arc. 

The linear-programming relaxation of MXFIR has integer extreme points, so we 
can take the dual of that relaxation to obtain MXFI-LRD after some simplifications 
(Derbes 1997): 


MXFI-LRD (A,R): 

z{A, R) = max - AR 

y 

s-t. X-f-t ■ Z-f-t =0 V/ G A 

k&FS(i) k^RS(i) 
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0 < } \/k € A 


Of course, z(l,R) is only a lower bound on z(R), but we can use this lower bound, 
and interdiction plans obtained during its solution to solve MXFI(i?), at least 
approximately. Proeedures to do this are presented in the next chapter. 
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III. NETWORK INTERDICTION BY LAGRANGIAN RELAXATION 

AND BRANCH-AND-BOUND 


This chapter explains our algorithm for solving MXFI. The algorithm uses MXFI- 
LRD as Bingol (2001) does, but adds two new procedures to decrease optimality gaps. The 
first procedure finds optimal Lagrangian multipliers, which Bingol’s approach does not do. 
The seeond procedure implements a restricted version of branch-and-bound. “LRD” and 
“MXFI-LRD” are used interehangeably below. 

A. PRELIMINARIES 

Derbes (1997) shows that Lagrangian relaxation can fail to solve MXFI(R) when 
capaeities are not unique, but may solve the problem eorrectly if the capacities are 
perturbed slightly so as to be unique. We use the proeedure suggested by Bingol (2001) to 
guarantee unique eapaeities: 

Procedure Perturb(u); 

{ 

Order u^. so that Wj < U 2 < ”• < u\a\ ; 

For (k and m such that ^k-l^ ^k+l ^k+m^ ^k+ni+l){ 

10"^ Uk, 

<—10 Mi + 1; 


10"^ Uk + m-1; 

} 

Return u; 


So, suppose that u = (1, 1, 3, 3, 3). After running this procedure, the network has 
unique arc capacities, u= (10000, 10001, 30000, 30001, 30002). These small perturbations 
are unlikely to change the optimal solutions to MXFI. (A scaling factor larger than lO"^ 
would have to be used if the network were very large and there were many arcs with 
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identical original capacities.) The use of this procedure is assumed and not shown 
explicitly below. 


B, THE LAGRANGIAN-RELAXATION MODEL 

The function z{X,R), from MXFI-LRD(/l,i?), is a piecewise-linear concave 
function in X for fixed R. Since LRD {X, R) relaxes MXFI-IP(i?), z(/l, R) < z{R) for all 
X , where z{R) is the optimal solution for MXFI-IP(i?). Naturally, we wish to maximize 
z{X,R) over/I. 

Given a fixed X , we can solve, or attempt to solve, MXFI(i?) and obtain the 
interdiction set {X) , for an unspecified R, through LRD {X, R) as follows (Bingol 2001): 

1. Modify G = (N, A) such that arc capacities become ul = mm{u^ > ^ A, 

and use a max-flow algorithm to find a minimum-capacity cut A ^, 

2. Define the interdiction set Ai{X) = {k g A^ \ u,^> /Ir^ }. 

The solution found through LRD(/l,R) identifies a minimum cut (A^,A,)that 
corresponds to a feasible or infeasible solution x to MXFI(R), where = 1 VA: g A^{X) , 
and = 0 otherwise. If x satisifies the budget constraint ^i;x^ < R , then it is a feasible 

keA 

solution. By solving LRD (/I, R) for several values of X, we may be able to find x, 
with^r^x^. = R . (Intuitively, this seems more likely to occur at X = u,^ which is why 

keA 

Bingol only evaluates such values of X.) In this case, MXFl(R) is solved exactly. If not, 
we will establish a non-zero optimality gap. 

To determine optimality gaps, we need upper and lower bounds on z(R). For a 
feasible interdiction set Aj{X) , the corresponding maximum flow is an upper bound (UB); 
z(/l,R)from LRD is a lower bound (LB) for z(R) for any/I; the best LB from the 
Lagrangian relaxation is found by solving max z{X,R) . 

A 

Derbes (1997) uses a binary search on the interval of uncertainty forX to maximize 
z{X,R) . His algorithm, when specialized for our version of MXFl(R), is: 
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Algorithm D(i?): 


Input: 


G= (N, A) directed network with source seN and sink teN; arcs indexed by k, 
u^. unique, scaled-up integer arc capacities > 0 , Vk e ^4, 


rk 

R 

Output: LB 


rk = \ for all interdictable arcs k e A; = <» otherwise, 

total interdiction resource available. 

lower bound for the network interdiction problem given R; 


{ 


} 


Initialize: <-oo ; A, <—m +l:/l <—m —l:/l<—/I : 

> 5 '^max max ’ mm mm ^ ^ '''max ’ 

while ) { 

u[ <— min|M^,/lr^} for all k ^ A; /* adjust arc capacities */ 

/* maxFlow( ) below returns a min cut A^ and max flow value / */ 
{Ac,f) <r- maxFlow (G, s, t, u'); 

Aj <— {k ^ Ac \ Uj^= ^ 

LB <— max {LB, f - XR) ; 

if (\Ai I = i?) { break; } /* from while loop */ 

lf{\A,\<R) 

if(l4l>i?) 


} 

print (“ LB for R =”, R, “is”, LB)-, 


In Algorithm D(i?), the interval of uncertainty initially has left endpoint 
4in = -1 and right endpoint +1. For ;L = + 4.„) / 2 , a max-flow 

problem is solved in G = (N,A) with u[ = max{u^,Xri^} and Aj is determined. The 
endpoints of the interval are adjusted by determining the sign of the slope of z{X,R) which 
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is sgn(^r^x^ \Aj\< R , A becomes the right endpoint, of the interval. If 

k 

\Ai\>R, then A becomes the left endpoint, A ^:^^, of the interval. If \Ai\ = R, then 
MXFI(i?) has been solved optimally. If an optimal solution to MXFI(i?) is not found, this 
process goes on until A is found such that A^:^^ <A< A^^ and A^^^ -A^-^ < 1. This A is 

close enough to X, because the arc capacities have been scaled up by a large 
multiplicative factor. (See Section A of this chapter.) 

With this procedure we can obtain a solution for only a single value of R. But, we 
are interested in solving MXFI(i?) for all “reasonable” values of R, which means 
|], where yfmin is a minimum-cardinality cut consisting of only interdictable 

arcs. (We assume such a cut exists.) In other words, \Amm\ is the minimum resource level 
that guarantees a feasible interdiction plan can force the network flow to zero. 

Performing a binary search for each i? e [0,| |] seems inefficient. This is 

especially true because z(/l,i?)and z{A,R') are closely related: 

z(A,R) = z(A,R') + A(R'-R). 

Furthermore, z(A,R) = z(A,0)-AR where z(A,0) is simply the maximum flow f{A) in 
G = {N,A) with arc capacities u[ = max{uj^,Ar^} \lk ^ A. 

Bingol (2001) attempts to solve MXFI(i?) for all (reasonable) values of R by 
approximately solving LRD(/l,i?) . He considers A in the range [0 ,m^^] and assumes the 

slope of z{A,R) = f{A)-AR can change only at A = Ui^. That is, he assumes that Aj{A) 
changes only in going from A = Ui^-£ to A = u,^+£ for certain and some small £■ > 0 . 
With this assumption, he can solve many instances of MXFI(i?) with acceptable accuracy 
by only evaluating z(/l,i?) at such points. Bingol’s algorithm is given next. The routine 
ProblematicRnO initially represents Bingol’s procedure ProblematicRl() for handling 
problematic i?-values. That will be replaced subsequently with our routine, 
ProblematicR2(), to improve results. 
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Algorithm A; Attempts to solve MXFI(i?) for eaeh value ofR g[ 0,| |] by computing 

z{X,R) for /I = and determining a feasible interdiction set 


Input; G={N,A) directed network with source A and sink teA; arcs indexed by A: 


rk 

Output: x{R) 

z{R) 

UB{R) 

LB(R) 


unique, scaled up integer arc capacities > 0 , \/k e A, 

Tyt = 1 for all interdictable arcs k e A; =co otherwise, 
optimal or suboptimal solution for MXFI(i?), 
interdiction set V i? e [0, |Amin|] that corresponds to x(i?) , 

for non-problematic R: optimal objective function value to MXFI(i?), 
for problematic R: feasible objective value, 
for problematic R: lower bound on optimal objective value, 
absGap{R) for problematic R\ absolute optimality gap UB{R) - LB{R), 
relGap{R) for problematic R\ relative optimality gap \QQxabsGap{R)ILB{R). 


/* define and order the candidate set of A, values */ 

Reorder arcs so that Wj < Wj < • • • < M|^| ; 

Initialize: R<— \A\+\; R' <— | ^ |+1; /<—0; /'<—0; /<—1; 

while (/ < \A \ and Rj^0){ 

A ^— Uj j 

u[ <— for all k e ; /* adjust arc capacities *! 

(Ac,/) <— maxFlow (G, 5, t, u'); 

Aj <— {k ^ Ac \ A}', 

4<—1VA:g4’ \/k € A\Aj; 

i? <— I y4/|; 
z <— /- AR; 

/* ProblematicRl() is called by Bingol; we will call ProblematicR2() */ 

if (R' -R > 1 ){ ProblematicRn(R',R,/',/,/l',/l,4); } 

if (R ^ R') { 
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print (“ The optimal solution to MXFI(i?) for R = ”, R, “is ”); 
print (x,Aj,z); 

} 

/* In the following ease, we know that MXFI(i?-l) will be solved by 
evaluating z(A + £,R-1) for a small s > 0 */ 
if ( u^. = X for some k' e ^/){ 

Ai<— Aj— {k'}; 

R^R-\; 

z z + ; 

Xj. <— 0; 

print (“ The optimal solution to MXFI(i?) for R = ”, R, “is ”); 
print (x,Aj,z); 

} 

if(R'^R) { 

R'<-R; A', <- A, ; z' <- z; 

} 

/<-/ + !; 

} 

} 


Algorithm A orders ares by increasing capacity Wj < Wj < • • • < M|^| and starts with 
/I = Mj. Because all arc capacities are A, Aj = A^, which corresponds to a post¬ 
interdiction flow of zero. This effectively defines the solution x(rj) to MXFI(rj) , where 
is the number of arcs in a minimum-cardinality cut consisting of only interdictable arcs. If 
= A for some k e A/, we also find a solution to MXFI(rj -1). If not, the next iteration 

sets A to and attempts to find a solution x(r 2 ) to MXFI(r 2 ) such that If 

r 2 = Tj -1, the algorithm continues without need of ProblematicRn(), and repeats for 
increasing values of A = u^. 
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When we eompute i? = | ^4^ , we know that we solve max z(A,IAj |) = z(/l*,i?), 

A 

i.e., we maximize the Lagrangian funetion z{A,R) for R = \Ai\. Furthermore, sinee 
=R, z(AI,R) = z(R); therefore we have solved MXFI(i?). However, diffieulties 

k^A 

arise when the algorithm finds z{Al,R) = z{R) and z{X*j^.,R') = z{R') for R'>R + \, but 
not for R + \,R + 2,...,R'-\. We eall these “missing” or “problematie i?-values.” For 
problematie i?-values, Bingol (2001) derives a heuristie, ProblematieRl() to find a feasible 
solution and lower bound LB on z(/l*,i?) . The feasible solution gives an upper bound UB 
for z{R) , so he ean eompute an absolute, or relative optimality gap whieh are defined as 
absGap=UB-LB and relGap= 100 x absGap/LB, respeetively. 


Procedure ProblematieRl ( R', R, /', /, A', A, ) 

Input: R eurrent R-value eorreetly solved, 

R' previous R-value eorreetly solved {R' > R + 1), 

A value of parameter that maximized z{A,R) to solve MXFI(R), 

A' value of parameter that maximized z{A,R') to solve MXFI(R'), 

/ maximum fiow^^), 

/' maximum flow / (A') , 

A'j interdietion set for MXFI( R'). 


Output: x(R) 

UB{R) 

LB{R) 

absGap(R) 

relGap{R) 


feasible (possibly optimal) solution for MXFI(R), 
interdietion set that eorresponds to x(R), 
feasible objeetive value for \(R ), 
lower bound on optimal objeetive value z{R), 
absolute optimality gap UB{R) - LB{R) for x{R) , 
relative optimality gap \00xabsGap{R)/LB{R) ) for x(R) . 


Initialize; UB - AR'; A^ Ap, 
for ( r = R' - 1 down to R + 1 ) { 
kmin<- argminn^; 

LB <—max{/- Ar,f'- A'r }; 

UB <— UB + Uk ; 

'^min 

absGap <— UB - LB', 
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relGap <— IQQxabsGapILB', 

A ^ ^min ’ 

Xj, <— 1 VA: e <— 0 \/k e AXAp, 

print (“ For problematic R = ”, r, “the approximate soln. to MXFI(i?) is: ”); 
print (x,Aj ,UB, LB, absGap, relGap ); 

} 

} 


Algorithm A requires the solution of a sequence of maximum-flow problems. We 
use the shortest augmenting path max-flow procedure, maxFlow(), (Edmonds and Karp 
1972), which finds the maximum flow/ and a minimum-capacity cut Ac- Since the flow 
on each arc is non-decreasing, this procedure uses the previous flow level from the last 
iteration, instead of starting from scratch. This is important for fast run-times. 

C. FINDING AN OPTIMAL FOR PROBLEMATIC R 

The procedure ProblematicRl() above does not find the best LB, i.e., it does not 
maximize z{X,R) for problematic 7?-values and therefore does not establish the best 
achievable optimality gap. We can attempt to find the best lower bound by exploring the 
region by defining a finer grid for X . But this is computationally expensive and 

requires solving many max-flow problems. At this point. Proposition 1 expedites the 
process. 

Proposition 1 (Derbes 1997): The Lagrangian function z{X,R) has a slope of 
I Aj (X) I -R , where | Aj (X) \ is the number of arcs interdicted using X. 
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Figure 1. NETSample. Letters and numbers represent are labels and eapaeities, respeetively. 
Unlabeled ares eannot be interdieted. All labeled ares have 4=1- 


Figure 1 shows a sample network, NETSample, whieh we will use to explain 
Proposition 1. We will eompute z(/l,l)and look at the funetion in detail; 



Figure 2. Lagrangian funetion for z(A, 1) for NETSample. The first number in parentheses is the 
value of parameter A and the seeond number represents the value of the funetion. 

If we eompute the slope msc of line BC in Figure 2, we see that ^ I I 
= 3-1=2. Likewise, =IAj(A)l-B = 1-1 = 0, whieh also gives us the optimal 
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solution for i? = l. (When the slope is zero, \ Aj{X)\ = R = \, and we obtain an exaet 
solution for MXFI(l).) These results are in aeeordanee with Proposition 1. 

For NETSample, if we run Algorithm A, we run into a missing value of R. We find 
a solution for i? = 3 at /I, = 3 and R = \ at \=1, but not for R = 2 . The solutions for 

R = 3 and i? = 1 are Af^ = {a, b, c} and Af^ = {e}. Bingol’s heuristie finds the UB and 
LB for i? = 2, with the following steps: 

i. 4 ^’ = Af^-{kmin) = {a,b,c)-{c) = {a,b), 

ii. UB = z(3) + Umin =1 + 7 = 8 , 

in. LB = max{/(/Ij) -2/lj,/(/I 3 )- 24 } = max {19 - 14, 10-6 } = 5, 

iv. absGap = 8-5 = 3, 

V. relGap =100 x 3/5 = 60%. 

The relative gap, relGap, is quite high. To reduee this optimality gap, we exploit 
Proposition 1 and traee the funetion z(A,R) more elosely. 
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Figure 3. Lagrangian Function of 7? = 2 for NETSample. The solid line denotes that the value of the 
function is computed for those d. values. The dotted line shows how procedure 
solveExact() computes the optimal T . 


ProblematicRl() finds LB = 5 for i? = 2, which is the value of the funetion at point 
E. But with Proposition 1, we ean look for over the region [3, 7] to improve LB and 
possibly solve MXFI. If y4/(/l^) = 2 (it is not in this ease), then MXFI(2) is solved 
exaetly. Otherwise, we obtain an improved LB. To aeeomplish this, we do the following: 

i. /* Find the slopes of assoeiated lines */ 

nicD - I (3) I - R = 3 - 2 = I, rrij^^ = | ^ 4 ^ (7) | -R = 1-2 = - 1, 

ii. /* ^2 is the optimal parameter value for R =2 */ 

z( X^,2) =4 + (1)( 3), z{X^,2)= 5+ {-!){X 2 - 1), 

iii. /* Compute from the equations above */ 

Xl = 5.5 and z{X*2,2) = 6.5, 

iv. /* Cheek for optimality */ 

AX5.5) j^2 , /* MXFI(2) is not solved exaetly */ 
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V. /* Compute the optimality gap with the new LB */ 
absGap = UB-z(X*2,2) = 8 - 6.5 = 1.5, and 
relGap = lOOxabsGap Iz{X^,2) = 100 x 1.5 / 6.5 = 23.1%. 


The relative gap has deereased from 60% to 23.1%. The best Lagrangian multiplier 
is ^2=5.5. Beeause Algorithm A eomputes z{X,R) at Ji = Ui^\/k e A, it would never 
diseover this LB. This example has only one missing i?-value. However, there ean be 
eontiguous problematie i?-values. The following proeedure will also handle these eases. 


Procedure ProblematieR2( R', R, f ,f, X', X,A'): 


Input: R 

R' 

A 

X' 

f 

r 


eurrent i?-value eorreetly solved , 
previous i?-value eorreetly solved (i?' > i? + 1), 
value of parameter at whieh R is solved, 
value of parameter at whieh R' is solved, 
maximum flow^ X ), 
maximum flow^ A! ). 


Output: x{R) 

4 ") 


optimal or suboptimal solution for MXFI(i?), 
interdietion set V i? e [0, |Amin|] that eorresponds to x(i?) , 


z(R) for non-problematie R: optimal objeetive funetion value to MXFI(i?), 

UB{R) for problematie R: feasible objeetive value, 

LB{R) for problematie R\ lower bound on optimal objeetive value, 
absGap{R) for problematie R\ absolute optimality gap UB{R) - LB{R), 
relGap{R) for problematie R: relative optimality gap \00xabsGap{R)/LB{R). 

For (r =i?' - 1 to R+\) { 

point <- X; prevPoint <- A'; R^^^^ <—R' ; <—R'; 


if(r = R + l){R_^R;} 
while {point prevPoint) { 
prevPoint <— point, 

point <— {f - f -XR + X'R')/{R' -R); P see the explanation below */ 
u[ <— minju^ , (point) } for all A: e 
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{Ac, flow) <— maxFlow(G, s, t, u'); 

AI <— {k ^ A^ I Uj^ — A,)', 

<^l\/k e Aj-, Xj.<—0 yk e A\Aj; 

/* If the exact solution is found, then it prints the results */ 
if (I 4 I = r) { 

print (“ The optimal solution to MXFI(i?) for R = ”, R, “is ”); 
print {x,Ai,z); 

prevPoint <— point ; /*get out of while loop */ 

if ( ^prev ~ Kur =1)1 ^prev ^ Kur ’ } 

} 

if (I 4 I > ^ ) { Al' <r- point; i?' <— | 4 I i f flow. A) <— Aj; } 
if(l4l<^){ point ; i? <— I 4 I i / ■^ flow,} 

} 

K ^ point ; /* The best Lagrangian multiplier is found */ 

if(4..v-4„.>i){ 

/* With revised values of X and A', we can reuse ProblematicRlQ to 
find a feasible solution and optimality gap for missing values of R */ 
ProblematicRl ( R^^^^, , /', /, A', A, ); 

/* Uncomment the next statement to invoke Branch() . Notice that UB{R) 
comes from ProblematicRl () */ 

!* Branch(4^^^, 44 -s,Af-^"\ UB{R)); */ 

D V_ D 

prev ^cur ’ 

} 

} 

} 

In our heuristic above, we derive the formula point = {f - f - AR + A' R') l{R' - R) 
in order to find the intersection with the following steps: 
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1. /* Find z(point,r) using the previous solution’s parameters (R',A',f'). 
Notiee that r is the missing value and slope is R' -r */ 

z{point, r) = (/' - /IV) + {R' - r){point - /I'), 

2. /* Find z(pomt,r) using the eurrent solution’s parameters (i?,/I,/). This 
time, slope is r-R */ 

z(point, r) = {f - Ar) + (r- R)(A - point ), 

3. /* Using these two equations, solve for point" */ 
point = {f-f'-AR + A'R') /(R' - R). 

When Algorithm A runs into one or more problematie i?-values, ProblematieR2() is 
ealled to find the best Lagrangian multiplier for eaeh missing value of R. This heuristie 
seeks over an interval of uneertainty by performing a binary seareh based on the sign of 

slope - R . First, it finds the interseetion, point, using the formula derived earlier. 

If {point) = r , then MXFI(r) is solved exaetly. Otherwise, it updates the endpoints as in 
Derbes’s proeedure (1997). It keeps updating until it eonverges to /I* and possibly finds 
an exaet solution. If there are multiple values of 7? e - i,Rcur +1] ’ ProblematioRl() is 

invoked to improve the results with the best Lagrangian multipliers for eaeh of these 
problematie 7?-values. For further improvements, we ean a invoke Braneh() whieh applies 
a restrieted braneh-and-bound proeedure. 

If there is only one missing 7?-value, one of two oases ooours, and it only takes one 
evaluation of z{A,R) to figure this out: 

1. If there is only a single breakpoint, then we must oompute it exaetly when 
we oompute the interseetion of the two lines extending upwards from the ooordinates 
{Ag,z{Aj^,,R')) and (/l^ ,z(/l^,7?)) in {A , z)-spaoe. Let that interseetion point be (A" ,z"). 

The single breakpoint is verified by oomputing z(A",R) and finding z” = z(A",R) . 

2. If there are two breakpoints in the interval {A^,,Ag), we oompute (A" ,z"), 
evaluate z(A",R) and find that z(A",R) < z". Consequently, the slope of the z{A,R) ourve 
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must lie strictly between one and negative one. Since the slope can only change by integer 
amounts, we have found a portion of the curve where the slope is zero and will, 
consequently, have solved MXFI(i?) optimally. 

If there are sequential problematic i?-values, there are three possibilities: 

1. There may not be an exact solution for any of the problematic i?-values. In 

this case, the procedure finds A* and then ProblematicRI() is called with updated 
parameters to find the best LB. In fact, A and A' computed earlier by Algorithm A are at 
some Uj^. However, when computed by ProblematicR2(), either A or X becomes the best 
Lagrangian multiplier, which helps improve LB for these problematic i?-values. To invoke 
Branch(), we use the following idea: When there are more than one problematic i?-values at 
one breakpoint, the best Lagrangian multiplier is the same for all these values. In other 
words, for i? e [i?' -1, i? +1], /l*^j = A^^^^ =••• = A]^_y. Therefore, we use only one of them 
with a slight perturbation, /l*^j - ^, to invoke Branch(), where | {A*_^^ - ^) | = = R'. 

2. The interval of uncertainty might yield an exact solution for some of the 
problematic i?-values. If it does, then the procedure prints out the results of those values, 
and Rp^^^ and establish new endpoints for the problematic i?-values. If R^^^^ -Rcur > 

ProblematicRIO and/or BranchQ can be invoked to find a feasible solution and optimality 
gap for problematic R e [R^^^ ^h^prev “1] explained in the previous paragraph. Since 

and/or might change for the problematic R e [R^^^ “1] ’ only LB 

but also UB are computed by ProblematicRl(). 

3. The procedure might solve MXFI(R) exactly for all missing R-values in the 
sequence. In our tests in Chapter IV, we observe this in two instances out of 10 test 
problems. 

D. RESTRICTED BRANCH-AND-BOUND 

Using procedure ProblematicR2() to solve MXFI, we may be able to reduce the 
optimality gap by a substantial amount or solve MXFI exactly, for certain values of R. 
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However, there might be some problematie i?-values for whieh the optimality gap is still 
large. For these i?-values, we ean resort to braneh-and-bound. 

The usual idea of the braneh-and-bound in the setting of integer programming is the 
following: Assume that we have a feasible, binary integer-programming problem (BIP). 
We ean solve the BIP as a linear program (LP) (i.e., the linear-programming relaxation of 
BIP) by ignoring the integer restrietions. If the LP optimal solution eontains only binary¬ 
valued variables we have solved the BIP. Otherwise, we seleet a fraetional variable, say 
x^, and braneh on it. That is, we ereate two new LPs, LPi and LP 2 , by introdueing the 

eonstraints = 0 and = 1, respeetively. An optimal solution to the original problem is 

eontained in the feasible region to at least one of these two problems. Next, we seleet one 
of the problems, say LPi, solve it (the LP relaxation that is), and evaluate its objeetive. If 
the objeetive is no better than the value of some ineumbent integer solution (whieh we 
obtained heuristieally, say), then that problem is fathomed by bounds. If the problem is 
infeasible, it is fathomed by infeasibility. If the solution is integer then we have fathomed 
by integrality (and we save this solution as the ineumbent if it is better than any we have 
seen before). If none of those eases oeeur for LPi, we seleet a fraetional variable and 
reeursively braneh on it. Braneh and bound stops when all nodes are fathomed: The best 
solution found is optimal. (Note that some subproblems that have not been solved may be 
fathomed using information from other subproblem solutions. For instanee, if an optimal 
integer solution to LPi is found and its objeetive value equals that of LP, then LP 2 is 
implieitly fathomed by bounding.) 

But in our seheme for solving MXFI, there are no fraetional values. Thus, we need 
a different partitioning seheme. We will base the partition on, essentially, setting = 0 or 

= 1 for A: e (/I') , where (A') > R. In this ease, “foreing are k into the solution” is 
equivalent to setting x^ = 1 and “foreing are k out of the solution” eorresponds to setting 
x^ = 0. And, we will eonsider a node fathomed if, at that node: 

1. \Af = R (equivalent to being fathomed by integrality; an optimal solution to 
the subproblem has been found), and 
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2. z(/l, R) > UB . (Fathomed by bounding.) 

If we use this teehnique in our problem, our enumeration tree might look like the 
following, given Aj {A,') = {m, n, r) : 



(Level 0) 


(Level 1) 


(Level 2) 


Figure 4. Sample Enumeration Tree for 7? = 1. Solutions are given next to nodes in curly braces and 
the branching decision is next to the branches: “ The initial solution at the root node at 
level 0 is {m, n, r). “m ” and “m” represent the exclusion and inclusion of arc m, 
respectively. Note that none of the nodes are fathomed in this sample enumeration tree. 


In Figure 4, every time we foree an are in, we obtain the same solution. For 
example, when are m is ineluded for the solution at Level 1, we get the solution {m, n, r) 
from Level 0. But, it is unneeessary to eonsider a braneh with {m, n, r}, beeause we have 
already evaluated this solution. So, we modify the proeedure as follows: 
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(Level 0) 


(Level 1) 
(Level 2) 


Figure 5. Modified Enumeration Tree for 7? = 1. Solutions are represented next to nodes in {}. mn 
represents the inelusion of arc m and exclusion of arc n. A full partition of the space of all 
possible solutions includes {m, n, r}, but this branch need not be created, because this 
solution has already been evaluated at the root level of enumeration tree. 

In this modified enumeration tree, we begin with the exclusion of the first arc m in 
the initial interdiction set and we obtain the solution {v,y} as a result. Since | ^4, | 1 for 

this solution, two new subproblems are introduced by adding the condition v or vy . This 
process continues for all nodes until each subproblem node is fathomed. When all the 
nodes are fathomed, the best feasible solution encountered is declared optimal. 

Because of limited time for this study, we restrict the enumeration tree to a depth of 
one (Level 1). So, for Aj = {k^,k 2 ,...}, LB is computed as: 

Z.5 = min{ max z(/l,i? I = 0), max z(/l,i? | =0),...}. 

A ^ A ^ ^ 

If the LB comes from a node at which \Aj\ = R, Aj is the optimal solution to MXFI(i?). 
This is true because all the other nodes are fathomed by bounding. 

We find feasible solutions for each branch created as follows: If \Ai\>R, we 
“uninterdict” arcs k e Aj, in order of increasing capacity, until \Ai\ = R. If | ^4, | < i?, we 
interdict arcs k e {A^, \Aj}, in order of decreasing capacity, until | y4J = i?. The feasible 
solution with the lowest LfB becomes the best feasible solution. 

If there is more than one missing value of R, we apply this procedure sequentially 
as follows ( A^ denotes the set of arcs to be included): 
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Procedure Branch( i?' , R, A, Aj^ ^, UB(r)) 


Input; 


Network G= (N, A), source node 5 , sink node t, 

R current i?-value correctly solved, 

R' previous i?-value correctly solved (i?' > i? + 1), 

A value of parameter such as A = A^^^ - s , 


A\^ ^ interdiction set that corresponds to x(i?'), 

UB{R) feasible objective value for \{R). 


Output: x{R) 

LB(R) 

UB(R) 

absGap(R) 

relGap{R) 


feasible (possibly optimal) solution for MXFI(i?), 
interdiction set that corresponds to x(i?) , 

possibly improved lower bound on optimal objective value z{R), 
possibly improved feasible objective value for x(i?), 
for problematic R: absolute optimality gap UB{R) - LB{R), 
for problematic R: relative optimality gap \00xabsGap{R)/LB{R). 


{ 

for (r = i?' - 1 to i? + 1) { /* Solve for all missing values of R */ 

u[ <— mm.{uj^,Aj^,rj^} for all k e ; /* adjust arc capacities */ 
LB{r) <— GO; 

a"" <— 0; /* This will be the set of arcs forced into the solution */ 
for (a = 1 to r + 1) { /* Include and exclude the arcs in AAr) */ 

Let k be the arc in AA r ); 

\i{k i A^) { /* if A: is not already forced in */ 

u[<r- Ui^; /* forces arc k out of the cut, i.e., sets x^-Q*l 

} 

{Ac,f) <— maxFlow (G, 5, t, u'); 

Aj {k e Ac I A}', 

A^ <- AjUA^; 

Xf. <—lVk e Ap Xf.<—0 VkeA\Ap 
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/* Find a feasible solution */ 
if(MJ>r){ 

for {\Aj\ down to r ) { 
^argmmn,; 


local ^ U^local + ^k„ 


A ^ A ~Amm} ’ 

<— 0; 


} 


} 

/* Find a feasible solution */ 
if(|^/ |<r) { 

for (r down to | ^4J ) { 


K 

UB 


<— arg max w. ; 

k^iAf-IA,} 


local ^ ^^local ’ 


4 <- ^/ + '■> 

<-l ; 

} 

} 

LBiocai =f-K'A-a + \)-, 


LB{r) = mm{LB{r),LB^^^J 
UB{r) = mm{UB{r),UB,^^J 

u[ = 0; /* Set x^ = 1 in the following subproblems */ 

A" <- A" vj{k)-. 


absGap (r) <— f/5 (r) - LB (r ); 
relGap (r) <— 1 OOxabsGap (r) ILB (r); 

print (“ For problematie i? = ”, r, “the approximate soln. to MXFI(i?) is: ”); 
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print (x,A^ ,UB (r), LB (r), absGap (r), relGap (r) ); 

} 

} 

In this procedure, we foree an are out of a solution by taking its eapaeity baek to its 
original level, > /l^,. This ensures that it will never be seen as interdieted sinee Aj = {k 
e Ac \ ul = . “Foreing in” is aeeomplished by setting u[=0, whieh means that are k is 

destroyed and eannot be used. 

Although it is not shown in the pseudo-eode for simplieity, \A^\^r, \nq keep 
solving max-flow problems for /I e to find | A^(A) \ = r. If it doesn’t exist, then we 
find and sueh that \Aj{A^.J\<r< \Aj{A^^^)\. After finding these values, we 
eompute the way it is eomputed in ProblematieRl( ): 

= max{ifiA^^) - (r - a +1)), (/(4„) - A ^^(r - a +1))} . 

If we LmALB = UB , then MXFI is solved exaetly for that value of R. Otherwise, we 
may obtain an improved LB and/or UB, whieh deereases the optimality gap. 
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IV. COMPUTATIONAL RESULTS 


In this chapter, we deseribe the design of test problems, and provide computational 
results for our Lagrangian-based proeedures to solve MXFI. 

A. TEST NETWORK DESIGN 

We use grid networks to test our procedures. “NET axb" denotes such a network, 
where a is the number of nodes on the horizontal axis and b is the number of nodes on the 
vertieal axis. The number of nodes in a grid network is \ N \ = ab + 2, and the number of 
arcs is \A\ = a{5b -A)-5b + 6 . A representation of a small network is given in Figure 6: 



Figure 6. NET 3x3. Arc capacities and resources are not shown in this network. The source 
node and sink node are numbered 1 and 2, respectively. 


As shown in Figure 6, all horizontal ares lie in the direction of “west to east,” and 
diagonal arcs are oriented in the southeast and northeast directions. There are also vertical 
arcs, north-pointing and south-pointing, between adjacent nodes in a eolumn of nodes; 
although these are omitted in the first and last eolumns, where they would be superfluous 
for maximizing flow. Ares going out of the source node and going into the sink node have 
infinite eapaeity and eannot be interdieted. (In fact, uninterdictable arcs have r^. = 999, 
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which effectively makes them uninterdictable.) All the other arcs are interdietable. Finite 
arc capacities are drawn from a diserete uniform distribution with range [1,50]. 

We use the network generator written in C by Derbes (1997) and modify it for 
different networks. The network files generated by this program provide the following 
information about the structure of the network: number of nodes, number of ares, source 
and sink nodes. Also, each arc k receives a capacity , cost of repair c^., resource needed 

to interdict , and “reparability” Pi^. However, we use only and in our tests. 

B, RESULTS 

We use six different grid networks for testing. The smallest one, NET 5x5, has 27 
nodes and 86 arcs. The largest one, NET 20x20, has 402 nodes and 1826 arcs. All tests 
are performed on a 533 MHz Pentium III personal computer with 192 MB of RAM and 
running the Mierosoft Windows XP Professional operating system. All programs are 
written in the Java 1.3 programming language. 

We use three different methods to test the effectiveness of our Eagrangian-based 
proeedures. Eaeh method uses Algorithm A with a different proeedure or procedures. 
Proeedure ProblematieRl() is the heuristie that Bingol (2001) uses to find a feasible 
solution, and a lower bound on z(Ag,R). We use ProblematieR2() to find the best 
Eagrangian multipliers. The proeedure Braneh() applies the restricted branch-and-bound 
technique introduced in Chapter III. Method I uses ProblematicRI() integrated into 
Algorithm A; this is the teehnique Bingol (2001) analyzes. In partieular, every time 
Algorithm A runs into a contiguous set of problematic 7?-values, ProblematicRl() is called 
to find UB and LB for eaeh. Method 2 works the same way, but ProblematicR2() is 
invoked to find for problematie 7?-values before parameters are passed on to 

ProblematicRlO which is, in this case, only used to find feasible solution. Eurthermore, 
the indieated eall to Branch() inside of ProblematioR2() is not made. Method 3 invokes 
proeedure ProblematieR2() exaetly as stated in the text, i.e., with a call to Branch() if a 
positive gap remains. 
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Since we are really only concerned with problematie i?-values, the tables only 
present: the problematie values R, total run-times, upper bounds z(i?), lower bounds 
z{X,R), absolute gaps “absGap,” and relative gaps “relGap.” Tables also inelude the 
uninterdieted maximum flow value /*. 

Table 1 shows results for Method 1. Although all absolute gaps are small, relative 
gaps can be as large as 25.0% when the divisor z{A,R) is small. Eaeh network that we 
have tested has yielded one or more problematie values. Run-times are small, the largest 
one being 31.5 epu seconds. 

Table 2 exhibits results for Method 2 using the same networks. Sinee the best 
Lagrangian mulitiplier is not necessarily an integer, z(/l*,i?) might have fractional values. 
“absGap” and “relGap” are the results computed with these fractional values. However, 
optimal objeetive values must be integral, so it makes sense to round z(A^,R) up to the 

nearest integer. Thus, “[absGap]” and “[relGap]” denote absolute and relative gaps 

eomputed after this rounding. Of 36 problematic i?-values found with Method 1, Method 2 
solves 16 exaetly. The largest remaining relative gap is 14.3% for R = 51 in NET 20x 20. 
Run-times increase only slightly for the networks that have a few problematic R-values. If 
ProblematicR2() is invoked many times, the inerease is about 25%. 

Table 3 presents the results of Method 3 which invokes proeedure Branch(). We 
only pass on to this proeedure the problematie R-values remaining after Method 2. With 
this method, the number of remaining problematic R-values values goes down to four for 
NET 20x20. The largest relative gap drops to 7.1%, sinee a better feasible solution is 
found by Branch(). There ean be, unfortunately, a drastie inerease in run-time; for 
example, the run time increases from 31.5 seconds to 4013.7 seconds for NET 20 x 20. 

To explore the effeetiveness of these methods further, we use NET 8x15 with 
restricted capacities. There are four types of restricted capacities, denoted 16[1,2], 8[1,4], 
4[1,8], 2[1,16]: The notation “/^[l,^]” indicates that the random eapaeities are drawn 
uniformly from the set {p,2p,...,qp}. Bingol (2001) finds problem instances like these, 
with few distinct capacities, to be particularly difficult. That is, they tend to yield large 
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optimality gaps. We apply the same three methods as before: Table 4 ineludes only 
problematie i?-values with relative gaps higher than 5% remaining after Method 1. 

We start with a total of 16 problematic values from Method 1. The largest relative 
gap from Method 1 is 20%. By using Method 2, the number of problematic i?-values goes 
down to eight, with the largest gap being 6.1%. Four of these problematic values are 
solved exactly by Method 3. However, there is no improvement in the relative gaps of the 
remaining four. 
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G ^(N, A) 

r 

R 

Run Time 
(cpu sec.) 

Problematic 

R -Values 

z(7?) 

z(A,R) 

absGap, 

relGap, 





6 

55 

53 

2 

3.8% 

NET 5x5 

234 

[1,13] 

0.1 










2 

159 

152 

7 

4.6% 





17 

18 

17 

1 

5.9% 





16 

29 

27 

2 

7.4% 





15 

39 

37 

2 

5.4% 

NET 8x 8 

416 

[1,22] 

0.4 










14 

48 

47 

1 

2.1% 





6 

202 

200 

2 






5 

228 

227 

1 


NET 8x 12 

690 

[1,34] 

1.0 

15 

189 

188 

1 

0.5% 





40 

5 

4 

1 






39 

7 

6 

1 

16.7% 





34 

27 

26 

1 

3.8% 





24 

146 

145 

1 


NET 8x 15 

841 

[1,43] 

1.9 










23 

165 

162 

3 

1.9% 





22 

181 

180 

1 

0.6% 





17 

279 

278 

1 

0.4% 





6 

590 

588 

2 

0.3% 





38 

11 

10 

1 






10 

431 

430 

1 

0.2% 

NET 15x 15 

878 

[1,43] 

7.5 










9 

463 

462 

1 

0.2% 





3 

686 

684 

2 

0.3% 





53 

8 

7 

1 

14.3% 





51 

16 

14 

2 

14.3% 






19 

18 

1 

5.6% 





41 

83 

80 

3 

3.8% 





40 

93 

89 

4 

4.5% 





39 

101 

98 

3 

3.1% 





38 

no 

108 

2 

1.9% 

NET 20x 20 

1173 

[1,58] 

31.5 

37 

119 

118 

1 

0.8% 





32 

175 

174 

1 

0.6% 





17 

508 

507 

1 

0.2% 





16 

542 

541 

1 

0.2% 





15 

576 

575 

1 

0.2% 





14 

611 

609 

2 

0.3% 





11 

712 

711 

1 

0.1% 





7 

859 

858 

1 

0.1% 


Table 1. Test Results for Method 1. This is essentially Bingol’s proeedure (Bingol 2001) whieh uses 

Algorithm A with ProblematieRl(), but not ProblematieR2() or BranehQ. f* shows 
maximum flow without interdietion. The eolumn labeled “R” gives the range of 
interdietion resouree levels tested; he maximum flow goes to 0 when R equals the seeond of 
the two numbers. Run times eover the whole proeedure, but only problematie R-values 
remaining after the proeedure is run are shown. z{R) is the upper bound obtained from the 
best solution with ProblematieRl() and z{X,R) is the largest value of Lagrangian funetion 
obtained from MXFI-LRD. z{A,R) is always integer sinee ^ = for ksA. 
“ absGapj ” is z(R) - z(X, R) and “ relGapj ” is 100 x absGapj / z(X, R). 
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G={N, A) 

/* 

R 


Problematic 

R -Values 

z{R) 

z(J.‘,R) 

absGap 

relGap 


absGap^ 

relGapj 

relGap, 





6 



1.5 

2.8% 

54 

1 

1.9% 

mmm 

NET 5x 5 

234 

[1,13] 

0.1 


■HH 

■ImH 






■mpIB 





2 


■Ml 

6.5 

4.3% 

153 

6 

3.9% 






17 

18 

17.2 

0.8 

4.7% 

18 

- 

- 

5.9% 





16 

29 

27.4 

1.6 

5.8% 

28 

1 

3.6% 

7.4% 





15 

39 

37.6 

1.4 

3.7% 

38 

1 

2.6% 

5.4% 

NET 8x 8 

416 

[1,22] 

0.7 














14 

48 

47.8 

0.2 


48 

- 

- 

2.1% 





6 

202 

200 

2.0 


200 

2 

1.0% 






5 

228 

111 

1.0 

■nKB 

111 

1 

0.4% 


NET 8x12 


[1,34] 

1.0 

15 

189 

188.5 

0.5 

0.3% 

189 

- 

- 

0.5% 





40 

5 

4.3 

0.7 

16.3% 

5 


- 

25.0% 





39 

7 

6.7 

0.3 

4.5% 

7 


- 

16.7% 





34 

27 

26.5 

0.5 

1.9% 

27 


- 

3.8% 





24 

146 

145.5 

0.5 


146 


_ 

0.7% 

NET 8x15 

841 

[1,43] 

2.3 














23 

165 

163 

2.0 

1.2% 

163 

2 

1.2% 

1.9% 





22 

181 

180.5 

0.5 


181 


- 

0.6% 





17 

279 

278.5 

0.5 


279 


- 

0.4% 





6 

590 

588.5 

1.5 


589 


0.2% 

0.3% 





38 

11 

10.5 



11 


- 

inn^H 





10 

431 

430.3 



431 


_ 


NET 15X 15 

878 

[1,43] 

10.3 














9 

463 

462.7 


0.1% 

463 


- 

0.2% 





3 

686 

684.5 

1.5 

0.2% 

685 


0.1% 

0.3% 





53 

7 

7 

- 

- 

7 

- 

- 

14.3% 





51 

16 

14 

2.0 

14.3% 

14 

2 

14.3% 

14.3% 





50 

19 

18 

1.0 

5.6% 

18 

1 

5.6% 

5.6% 





41 

83 

80.5 

2.5 

3.1% 

81 

2 

2.5% 

3.8% 





40 

93 

90 

3.0 

3.3% 

90 

3 

3.3% 

4.5% 





39 

101 

99.5 

1.5 

1.5% 

100 

1 

1.0% 

3.1% 





38 

110 

109 

1.0 

0.9% 

109 

1 

0.9% 

1.9% 

NET 20X20 

1173 

[1,58] 

39.3 

37 

119 

118.5 

0.5 

0.4% 

119 


- 

0.8% 





32 

175 

174.5 

0.5 

0.3% 

175 


- 

0.6% 





17 

508 

507.8 

0.2 


508 


- 

0.2% 





16 

542 

541.5 

0.5 

0.1% 

542 


- 

0.2% 





15 

576 

575.3 

0.7 

0.1% 

576 


- 

0.2% 





14 

609 

609 

- 

- 

609 


- 

0.3% 





11 

712 

711.5 

0.5 

0.1% 

712 


- 

0.1% 





7 

859 

858 

1.0 

0.1% 

858 


0.1% 

0.1% 


Table 2. 


Test Results for Method 2. This method uses Algorithm A and invokes ProblematieR2() 

* 

and ProblematieRl(). ProblematieR2() finds the best Lagrangian multipliers for 
problematie R-values, whieh are then passed on to ProblematieRl() to find feasible 


solutions. Sinee z{R) is integer, 


<\,R) 


is a valid lower bound on z{R). “ absGap 2 ” is 


z{K)- 




and “relGap 2 ” is 100 xabsGap 2 / 




‘ absGap' 


IS 


z(R)-z(X ,R) and “relGap” is 100 x absGap / z(/l ,R). indieates that MXFI is 

R R 

solved optimally. Only two problematie R-values, R = 53 and R = \A for NET 20 x 20 , are 

* 

solved exaetly before rounding up z{l„,R) but, overall, more than 50% of the problematie 

R 

R-values are solved with rounding. relGap^ is shown for eomparison. 
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G =(N, A) 

* 

R 

Run Time 

Problematic 



\LB^ 

absGapj 

relGapj 

relGap, 

relGapj 

/ 

(cpu sec.) 

R -Values 

z(R) 

LB 

NET 5X5 

234 

[1,13] 

0.2 

6 

55 

55 

55 

- 

- 

3.8% 

1.9% 

2 

159 

158 

158 

1 

0.6% 

4.6% 

3.9% 





17 

18 

18 

18 

- 

- 

5.9% 

- 





16 

29 

29 

29 

- 

- 

7.4% 

3.6% 





15 

39 

39 

39 

- 

. 

5.4% 

2.6% 

NET 8x8 

416 

[1,22] 

7.0 

14 

48 

48 

48 

- 

- 

2.1% 

- 





6 

202 

202 

202 

- 

- 

1.0% 

1.0% 





5 

228 

228 

228 

- 

- 

0.4% 

0.4% 

NET 8x12 

690 

[1,34] 

4.4 

15 

189 

189 

189 

- 

- 

0.5% 

- 





40 

5 

4.3 

5 

- 

- 

25.0% 

- 





39 

7 

6.7 

7 

- 

- 

16.7% 

- 





34 

27 

27 

27 

- 

- 

3.8% 

- 

NET 8x15 

841 

[1,43] 

25.2 

24 

146 

146 

146 

- 

- 

0.7% 

- 



23 

165 

165 

165 

- 

- 

1.9% 

1.2% 





22 

181 

181 

181 

- 

- 

0.6% 

- 





17 

279 

279 

279 

- 

- 

0.4% 

- 





6 

590 

590 

590 

- 

- 

0.3% 

0.2% 





38 

11 

11 

11 

- 

- 

MM 

- 

NET 15x15 

878 

[1,43] 

96.0 

10 

431 

431 

431 

- 

- 

Wm 

- 




9 

463 

463 

463 

- 

- 

mm 

- 





3 

686 

686 

686 

- 

- 

0.3% 

0.1% 





51 

15 

14 

14 

1 


14.3% 

14.3% 





50 

19 

18 

18 

1 

5.6% 

5.6% 

5.6% 





41 

83 

83 

83 

- 

- 

3.8% 

2.5% 





40 

93 

92 

92 

1 

1.1% 

4.5% 

3.3% 





39 

101 

101 

101 

- 

- 

3.1% 

1.0% 





38 

no 

no 

no 

- 

- 

1.9% 

0.9% 

NET 20x20 

1173 

[1,58] 

4013.7 

37 

119 

119 

119 

- 

- 

0.8% 

- 





32 

175 

174.5 

175 

- 

- 

0.6% 

- 





17 

508 

508 

508 

- 

- 

0.2% 

- 





16 

542 

542 

542 

- 

- 

0.2% 

- 





15 

576 

576 

576 

- 

- 

0.2% 

- 





11 

712 

712 

712 

- 

- 

0.1% 

- 





7 

859 

859 

859 

- 

- 

0.1% 

0.1% 


Tables. Test Results from Method 3. Method 3 is essentially Method 2 except that Branch() is 
invoked for problematic R-values that Method 2 cannot resolve. LB - min 
{max z(/l,R I = 0), max z(/l,R | “O)’---}- “absGap 3 ”is z(/l) - [~Z.R~] and 

“ relGap 3 ” is 100 x absGap 3 / f LB^. Absolute and relative gaps computed with fractional 

values are not shown in this table. 12 problematic R-values out of the 16 remaining after 
Method 2 are solved exactly by Method 3. UB for R = 51 of NET 20 x 20 drops from 16 

to 15, because a better feasible solution is found by Branch(). relGap^ and relGap 2 are 
also shown in the table for comparison. 
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Method 1 

Method 2 

Method 3 

GHN,A) 

Problematic 

R - Values 

z{R) 

z(A,R) 

relGapj 

z(R) 


relGapj 

z(R) 

\LB^ 

relGapj 


24 

320 

304 

5.3% 

304 

304 

- 

- 

- 

- 


23 

352 

320 

10.0% 

320 

320 

- 

- 

- 

- 


22 

384 

336 

14.3% 

336 

336 

- 

- 

- 

- 

16[1,2] 

21 

400 

352 

13.6% 

352 

352 

- 

- 

- 

- 

20 

416 

368 

13.0% 

368 

368 

- 

- 

- 

- 


19 

432 

384 

12.5% 

400 

394 

1.5% 

400 

394 

1.5% 


18 

448 

400 

12.0% 

432 

420 

2.9% 

432 

420 

2.9% 


17 

464 

432 

7.4% 

464 

446 

4.0% 

464 

446 

4.0% 


31 

104 

96 

8.3% 

96 

96 

- 

- 

- 

- 

8[1,4] 

30 

120 

104 

15.4% 

112 

no 

1.8% 

112 

112 

- 


29 

128 

120 

6.7% 

128 

123 

4.1% 

128 

123 

4.1% 

4[1,8] 

38 

240 

200 

20.0% 

200 

200 

- 

- 

- 

- 

37 

280 

240 

16.7% 

240 

240 

- 

- 

- 

- 


35 

280 

260 

7.7% 

280 

264 

6.1% 

280 

280 

- 

2[1,16] 

34 

340 

320 

6.3% 

340 

329 

3.3% 

340 

340 

- 


33 

400 

380 

5.3% 

400 

393 

1.8% 

400 

400 

- 


Table 4. Comparative Test Results for NET 8 x 15 with Restricted Arc Capacities. “Methods” are defined in the text. Arc capacities denoted p\l,q\ 

are drawn from the SQi { p,2p,...,qp }. “relGap^,relGap^, and relGap^ ” are the relative gaps obtained from Method 1, 2 and 3, respectively. 

Only R-values with relative gaps remaining greater than 5% after Method 1 are shown in this table. 50% of problematic R-values are solved 
exactly by Method 2, and 50% of those remaining are then solved by Method 3. Notice that UB for R = 30 of the network with arc capacities 
8[1,4] has dropped from 120 to 112, because R=31 is solved exactly by Method 2, and this makes it possible to find a better feasible solution 
for R = 30. 
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V. CONCLUSIONS AND FURTHER RESEARCH 


A. SUMMARY 

The maximum-flow network-interdiction problem (MXFI) arises when an 
“interdictor” uses limited interdiction resources to restrict an adversary’s use of a network. 
Conceptually, the interdictor destroys arcs to minimize the maximum flow that the 
adversary, or “network user,” can achieve through the network. In this study, we have 
explored new Lagrangian-relaxation techniques for solving MXFI for all values of 
interdiction resource i? in a specified range. We assume that one unit of resource is 
required to interdict an arc, so we only deal with integer values of R. 

MXFI can be solved by means of integer programming, but this approach is slow. 
Bingol (2001) modifies the Lagrangian-relaxation technique developed by Derbes (1997) 
to solve MXFI for all R e [l,|^,„i„ |], where denotes a minimum cardinality cut 

consisting of only interdictable arcs. The Lagrangian-relaxation technique solves a 
sequence of max-flow problems with related arc capacities, and is therefore quite efficient. 

Let z(R) denote the optimal objective value to MXFI for a given R, and let 
z{A,R) denote the Lagrangian lower bound on z(R) given Lagrangian multiplier A. 
Essentially, Bingol (2001) evaluates z(A,R) simultaneously for all R, but only for values 
of A equal to perturbed arc capacities. We follow this procedure—the range of A is 
correct—^but, when z{A,R) <z{R) for /l*= argmax z(/l,7?), where z{R) is the objective- 

A 

function value for the best feasible solution found. This is accomplished with a search 
technique based on the slope of z(A,R) with respect to A . For most test problems, we can 
decrease the optimality gap for “problematic 7?-values” by approximately 50%. A value of 
R is “problematic” if Bingol’s procedure leaves z{A,R)<z{R) . Although this approach 
focuses on tightening lower bounds, it is also finds better feasible solutions on occasion, 
and sometimes proves optimality of the best solution. However, some large optimality 
gaps may remain. 
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We tackle the remaining problematic i?-values with a modified branch-and-bound 
procedure. This procedure uses a partition based on the best interdiction set found, 
^4/= . In particular, we solve max z(/l,i? | = 1,=0),..., 

/I ^ ^ 

max z(/l, i? I x^_ = = k • • • r , = k = 0) . (max z{X, R\Xj^ = \\/k ^ A^) is already 

solved.) The minimum lower bound from these problems is a lower bound on the optimal 
solution. This partitioning scheme can be applied recursively to solve MXFI exactly for a 
given R, but we have only implemented a single application of the partition. That is, we 
have limited the enumeration tree of this method to one level below the root. Using this 
technique, we may find a lower bound which equals the upper bound obtained from the 
best feasible solution. In this case, MXFI is solved. In our tests, we have solved 
approximately 70% of the problematic 7?-values that remain after applying the first 
technique. But, this procedure runs slowly because of the need to solve many max-flow 
problems, so it is best used only when large optimality gaps are encountered. 

Having used these two techniques together, we have been able to solve MXFI either 
exactly or with small optimality gaps for most of the problematic 7?-values encountered. 
But, we still cannot guarantee an exact solution or a small optimality gap. 

B, FURTHER RESEARCH 

In this thesis, we have used a restricted modified branch-and-bound procedure that 
explores the enumeration tree only one level below the root. That is why some problematic 
7?-values still remain. A fully recursive implementation of the branch-and-bound algorithm 
will solve MXFI exactly, at least in theory. Even the restricted branch-and-bound 
procedure becomes slow because of the need to solve many max-flow problems. Thus, a 
fully recursive procedure based on this methodology would be even slower. To make such 
a procedure reasonably efficient, a faster maximum-flow algorithm should be used and 
techniques for reoptimizing related maximum-flow problems should be explored. 

For the sake of simplicity, we have assumed that that = 1 for all interdictable arcs 
k. That is, one unit of resource is needed to interdict any interdictable arc. Allowing 
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> 1 (but still integer) will make the problem harder. In this case, a complete branch-and- 
bound procedure would become more important to possess. 

We have assumed that grid networks are adequate representatives for real-world 
networks, so we have performed our tests on these networks. However, we are uncertain 
how well these procedures will perform on actual networks. Obviously, future research 
should apply these techniques to more realistic networks. 
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